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ABSTRACT 

Aims. We report simultaneous observations of the nearby flare star Proxima Centauri with VLT/UVES and XMM-Newton over three nights in 
March 2009. Our optical and X-ray observations cover the star's quiescent state, as well as its flaring activity and allow us to probe the stellar 
atmospheric conditions from the photosphere into the chromosphere, and then the corona during its different activity stages. 
Methods. Using the X-ray data, we investigate variations in coronal densities and abundances and infer loop properties for an intermediate- 
sized flare. The optical data are used to investigate the magnetic field and its possible variability, to construct an emission line list for the 
chromosphere, and use certain emission lines to construct physical models of Proxima Centauri's chromosphere. 

Results. We report the discovery of a weak optical forbidden Fe xm line at 3388 A during the more active states of Proxima Centauri. For the 
intermediate flare, we find two secondary flare events that may originate in neighbouring loops, and discuss the line asymmetries observed 
during this flare in Hi, Hei, and Can lines. The high time-resolution in the Ha line highlights strong temporal variations in the observed 
line asymmetries, which re-appear during a secondary flare event. We also present theoretical modelling with the stellar atmosphere code 
PHOENIX to construct flaring chromospheric models. 

Key words, stars: activity - stars: magnetic fields - stars: chromospheres - stars: coronae - stars: late-type - stars: individual: Proxima Centauri 
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1. Introduction 

Coronal heating of cool stars is thought to be driven by mag- 
netic fields generated in the stellar convection zones. For stars 
from spectral type F to mid-M, the magnetic fi eld productio n 
process can be modelled with an aQ, dynamo ( Par ked 1 1955b . 
How ever, stars of spectral type M3 or later are fully convec- 



tive (Chabrier& Baraffe 1997; Dorman et al. 1989) and are 



not expect ed to undergo the same dynamo process as earlier- 
type stars ( Browning] 2008 ). Other processes such as turbu- 
lent or a 2 dynamos should not be as effective as the aQ dy- 
namo and, indeed, many late M dwarfs ha ve rather low X-ray 
lumin osities in their states of quiescence (IRobrade & Sch mitt 



2005). Nevertheless, the same stars are capable of producing 



strong flares of short as we ll as long durations dStelzer et al 
2006t IRobrade etalj 1201 d) . so there must be some mech- 



Send offprint requests to: B. Fuhrmeister 

* Based on observations collected at the European Southern 
Observatory, Paranal, Chile, 082.D-0953A and on observations ob- 
tained with XMM-Newton, an ESA science mission with instruments 
and contributions directly funded by ESA Member states and NASA. 
** Full Table[6]is only available in electronic form at the CDS 



anism to allow violent releases of large amounts of mag- 
netic energy. Magnetic flux densities are enormous in mid- 
to late-M stars, providing firm evidence of an efficient mech- 
anism to produce and main tain strong magnetic fields (e.g. 
Reiners & Basil d2007l |2010)). In summary, the magnetic ac- 



tivity processes in late M dwarfs are not yet thoroughly un- 
derstood in detail. Multi-wavelength observations cover dif- 
ferent activity indicators and can therefore help us to recon- 
struct a picture of the whole stellar atmosphere. Examples of 



multi-w a velength obse r vation s c an be found in e . g.lOsten et al 



(200 5al). iBerger et all d2008l) . iKowalski et alj d2010al) . and 
iFuhrmeister et al.1 d2008l) . 

The energy released in flares is thought to originate from 
the interaction of magnetic fields in the convective zone of 
the star and photospheric motions that entangle the magnetic 
loops penetrating the stellar surface. The energy release is trig- 
gered by magnetic instabilities that cause the entangled mag- 
netic field lines to reconnect in the corona leading to heating 
and particle acceleration in the reconnection region. Particles 
are accelerated downwards into the chromosphere where they 
collisionally heat the denser plasma, which in turn expands and 
evaporates into the corona. The coronal loop is filled by dense, 
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hot plas ma, which is then ob servable as a flare in soft X-rays 
(see e. g. lHaisch et~ai1dl99ll) ). 

Proxima Centau ri, with a distance of only 1.3 pc 



(Ivan Leeu wen 2007) the star closest to the Sun, is a magnet- 
ically active star of spectral type dM5.5. It has be en frequently 
obser ved by various X-ray satellites: Ei nstein ([Haisch et alJ 
19801) . EXOSAT- IUE (IHaisch et all 1 19831) . XTE dHaisch et al. 
1998), ROSAT dVoges et al.1 1 1999) . A SCA dHaisch et alJ 
1995b. Chandra dWargelin & Drake! I2002T) . and XMM-Newton 
(Gu del et al.ll2004l) . Its quiescent X-ray luminosity varies in the 
range Lx ~ (4 - 16) x 10 26 erg s _1 , which is comparable to that 
of the Sun despite its 50 times smaller surface area. Over the 
past 30 years, several X-ray flares of Proxima Centauri have 
been observed, with the most extreme pe ak luminosities ob - 
served in a 2001 XMM-Newton observation dGiidel et alJ2004 . 
which exceeded typical quiescent state X-ray fluxes by a fac- 
tor of « 100. Here we present new data for Proxima Centauri's 
coronal and chromospheric properties, together with simulta- 
neous measurements of its large-scale magnetic field strength. 

Our paper is structured as follows. In Section 2, we describe 
our observations obtained with VLT/UVES and XMM-Newton. 
In Section 3, we compare the timing behaviour of Proxima 
Centauri in different wavelength bands. The coronal properties 
of Proxima Centauri such as temperatures and elemental abun- 
dances are presented in Section 4, while Section 5 describes 
the chromospheric and transition region properties of the star. 
Section 6 and 7 contain a discussion of the presented findings 
and our conclusions. 



2. Observations and data analysis 

The multi-wavelength observations reported in this paper were 
obtained strictly simultaneously with XMM-Newton and ESO's 
Kueyen telescope equipped with the Ultraviolet- Visual Echelle 
Spectrograph (UVES) on 9, 11, and 13 March 2009 (labelled 
"night 1", "night 2" and "night 3" in the following; see also 
Table[T]i. Due to the proximity of Proxima Centauri, interstellar 
absorption is negligible for the optical as well as for the X-ray 
data. 



2.1. Optical UVES data 

For our optical observations, the UVES spectrograph was op- 
erated in a dichroic mode leading to a spectral coverage from 
about 3290 A to 4500 A in the blue arm and 6400 A to 
10080 A in the red arm with a small gap from 8190 A to 
8400 A caused by the CCD mosaicQ. For the red arm, a non- 
standard setup was used to ensure coverage of the Ha line. We 
used exposure times varying from 1000 to 1800 seconds for 
the blue arm and of 90 to 450 seconds in the red arm due to 
variable seeing conditions. In the blue arm, we obtained 24, 
17, and 15 useful spectra during the three nights and in the red 
arm 215, 168, and 179 useful spectra, respectively. The typical 
resolution of our spectra is ~ 45 000. The red ar m spectra were 
reduced using the UVES pipeline vers. 4.3.0 (Ball ester et al 



20000. The blue arm spectra could not be properly reduced 
with the pipeline software but had to be r educed manually us- 
ing th e IDL reduction software REDUCE dPiskunov & Valenti 
20021) for several reasons. To obtain a simultaneous coverage 
with XMM-Newton, the star had to be observed at rather high 
airmass starting with a maximum of 2.4. Since the positioning 
of the star was accomplished with the red arm, a wider slit had 
to be used for the blue spectra, which resulted in an overlap of 
the bluest spectral orders. Therefore no dark exposure could be 
acquired and the flat-field and science spectra extraction had to 
be carried out with a fixed width and without scattered light. To 
complicate the spectral reduction even more, quite a number of 
spectra are contaminated by solar stray-light (from the Moon), 
which is outshone by the star in the red part of the spectrum but 
can be clearly recognised in the blue part of the spectrum. 

The wavelength calibration in the blue arm was carried 
out with thorium-argon spectra with an accuracy of ~ 0.04 A. 
Since the weather conditions and especially the seeing varied 
throughout the nights, no absolute flux calibration with a stan- 
dard star could be performed. To obtain an a posteriori flux cal- 
ibration, we used synthetic stellar spe ctra provided by the stel 



lar atmosphere program PHOENIX (Hauschildt et al. 1999). 



We determined the best-fit stellar model using a grid with T e g 
varying from 2700 K to 3400 K in steps of 100 K and log g of 
4.5, 5.0, and 5.5 and found that model to have T eff = 3100 K 
and log g - 5 . 5 in g ood agreement with the values derived by 
Dem orvetaf] (120091) . who had derived r eff = 3098 + 56 K and 



log g — 5.2 using the VLTI. For the blue spectra, only data at 
wavelength redward of 4000 A were used for the calibration 
since at short wavelengths the data are dominated by chromo- 
spheric emission and thus the model spectra are unreliable. We 
estimate that our errors in the flux calibration in the blue band 
are about a factor of two. A comparison wit h a flux-calibrated 



spectru m o f Proxima Centauri acquire d by ICincunegui et al 



(2007) and ICincunegui & Mauasl (2004) shows that our fluxes 
are higher by about a factor of two for the blue band and by a 
factor of two to three in the red band. Since these authors used 
a standard star and low resolution spectra for their flux calibra- 
tion, the main source of error in our calibration seems to be the 
model spectrum. 

In addition to the spectral data, we obtained a blue and a 
red light curve using the UVES exposure-meters, i. e., one blue 
and red photometer located in the two arms of the spectro- 
graph. These data are normally obtained for engineering pur- 
poses only and are not flux calibrated. We note that the flux 
bands of the exposure-meters are not identical to the spectral 
bands, the blue flux band in particular being "redder" than the 
blue spectrum. 

2.2. X-ray data 

Each of our VLT observations was accompanied by three si- 
multaneous 30 ks observations conducted with XMM-Newton; 
the exact observation times with the ObsIDs being given in 
Table [T] On board XMM-Newton, there are three telescopes 



1 A detailed description of the UVES spectrograph is available at 
http://www.eso.org/instruments/uves/doc/ 



2 The UVES pipeline manual can be found at 
ftp.eso.org/pub/dfs/pipelines/uves-cpl/uves-pipeline-manual- 1 3 .0.pdf 
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Table 1. The three XMM-Newton observations of Proxima 
Centauri in March 2009. 



ObsID 


start time 


duration (ks) 


label 


0551120301 


2009-03-10 02:23:34 


28.7 


night 1 


0551120201 


2009-03-12 02:11:19 


30.7 


night 2 


0551120401 


2009-03-14 02:20:45 


28.7 


night 3 



focusing X-rays onto three CCD cameras (one PN and two 
metal oxide semi-conductor (MOS) cameras with a sensitivity 
range of « 0.2 - 15 keV), which together form the European 
Photon Imaging Camera (EPIC). The X-ray telescopes with 
the MOS detectors are also equipped with reflection gratings. 
The two RGS (Reflection Grating Spectrometers) provide high- 
resolution X-ray spectroscopy (E/AE « 200-800) in the energy 
range of 0.35-2.5 keV capable of resolving individual X-ray 
emission lines. The X-ray instruments are accompanied by the 
Optical Monitor (OM), an optical/UV telescope that can be 
used with different filters for imaging and time-resolved pho- 
tometry. Useful data of Proxima Centauri were obtained with 
the OM, EPIC, and RGS detectors, which were all operated si- 
multaneously. The PN and MOS detectors were operated with 
the medium filter in full frame and large window mode, respec- 
tively. The OM was operated in fast mode with 0.5 s cadence 
using the U band filter covering a band pass of 300 - 390 nm. 

All XMM-Newton data were reduced using the standard 
XMM-Newton Science Analysis System (SAS) software, ver- 
sion 10.0. The EPIC light curves and spectra were obtained 
using standard filtering criteria. We adopted an extraction ra- 
dius of 15" centred on the source; for background subtraction 
we used nearby source-free areas. Source counts with energies 
in the energy range 0.2 - 10 keV were considered for the sci- 
entific analysis. For the first two nights, the X-ray observations 
displayed only a low background level. During the third night, 
short phases with high background levels were present; these 
time intervals were excluded from our spectral analysis. The 
X-ray light curves are background-subtracted and binned by 
100 s unless stated otherwise. 

Spectral an alysis was carried out with Xspec V12.5.0 
(Arnaud 1996 ) for the ov erall fitting processes and CORA 
dNess & Wichmann 2002) ? or fitting of individual spectral 
lines. For the overall analysis in Xspec, we used models 
with several temperature components assuming the same el- 
emental abundance for each component. These models are 
based on a collisiona lly ionised optical ly thin gas calculated 
with the APEC code dSmith et al.ll200ll) . Abundances are cal- 
culated relative to the s olar photospheric values given by 
Grevesse & Sauvalidl998h . 



3. Multi-wavelength timing behaviour 

In Figures [TJ [2] and [3] we present the background-subtracted 
X-ray light curves (taken with EPIC-PN) and the OM U-band 
light curve for our three Proxima Centauri observations, as well 
as the optical UVES blue-band flux and several chromospheric 
line fluxes; the X-ray light curves are binned in units of 100s, 



and data for the OM in units of 10s, while the UVES exposure- 
meter data are binned in units of 5s. 



3.1. Quiescence and small flares 

In the optical (UVES exposure-meter and in the OM), flar- 
ing state and quiescent phases can easily be identified even 
for smaller flares. For the chromospheric emission line light 
curves and the X-ray light curves this is not so easy: the lowest 
parts of the X-ray light curve - as seen in Figs. [TJ-[3]- show 
a typical count rate of « 1 cps in the PN detector, comparable 
to the low count rate states found in the 2001 observations of 



Proxima Centauri (see Giidel et all (|2002)). During these time 
intervals, the light curves also display some variability, i. e. do 
not remain constant. To be consistent with other authors, we 
nevertheless use the term "quiescence" and define it for our 
subsequent spectral analysis as those times with a PN count 
rate below 1.8 cps. 

In addition to the low-level variability seen during quies- 
cence, small flares are observed at 3:30 and 6:00 UT in the X- 
ray data from night 1 (see Fig. [TJ. Similarly, in the X-ray data 
from night 2 (Fig. O, which have a slightly higher mean count 
rate, there is some low-level variation with small flares at about 
4:40, 6:15 and 8:40 UT. Most of these flares are also visible in 
the OM light curves, while the UVES exposure-meter does not 
show significant variability. Apart from these small flares, the 
X-ray light curves also show some plateau-like structures and 
slow rises and decreases not seen in the OM light curve, which 
displays a constant quiescent level during these times. 

Nevertheless, a comparison of the X-ray flux with several 
chromospheric line fluxes (see Figs. [TJ and |2| shows quite good 
agreement, especially for the Ha light curve, which can be 
studied with the much higher time-resolution of the red arm 
spectra. For instance, the first Ha flare in night 1 at about 
3:30 UT can also be identified in the X-ray light curve, al- 
though not in the optical (UVES) light curve. In addition, 
the flare in night 2 at 6:15 UT is found in X-ray, OM, and 
in the chromospheric emission lines, but is much more pro- 
nounced in the Can H & K line than in the Balmer lines. This 
shows that the chromospheric emission lines can also be af- 
fected by small flares, when the continuum exhibits no changes. 
Moreover, the chromospheric emission lines exhibit variability 
that is not seen in the X-ray band, thus should be confined to 
the chromosphere; for instance, in night 2, during the first two 
hours, the Ha and other chromospheric emission lines display 
a flux decline, but the X-ray level is constant, and in night 1, 
the chromospheric lines show a peak at about 7:00 UT that is 
not noticed in the X-ray or c ontinuum. This beh aviour during 
flares has been found before. Osten et al. (2005b) found at ra- 
dio, optical, and X-ray wavelengths for the flare star EV Lac 
flares with no cou nterparts at other wavelengths. In addition 
Hilton et al. (2010) searched for flares in the SDSS and found 
243 flaring spectra of 63 flare stars with only two having an en- 
hanced continuum. Continuum en hancement occurs espec ially 
for the strongest flares (see e. g. iKowalski et alj (|2010a) and 
Hawlev&Fisnerl(ll992l) ). 
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Fig. 1. Light curves of Proxima Centauri from night 1 as seen 
in the fluxes of characteristic chromospheric emission lines 
(top, note the better time binning provided by the Ha line) as 
well as optical UVES, OM, and X-ray light curves (bottom). 
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Fig. 2. As Fig. [T] but for night 2. The two gaps in the UVES 
light curve are due to technical problems that led to the loss of 
the guide star. The legend is the same as in Fig.Q] 



We also analysed Proxima Centauri's integrated flux using 
PN spectra from the three individual XMM-Newton observa- 
tions. The X-ray flux levels did not change significantly for 
the first and second exposure, whereas there is a significant 
change in the X-ray flux in the third exposure as can be seen 
in Fig.|4l here the X-ray flux was determined using the Xspec 
best-fit spectral models (see Section 4) integrated in the energy 
range 0.2 to 10.0 keV. Specifically, the night 1 EPIC PN expo- 
sure results in an X-ray flux of 3.0 x 10 -12 ergcirT 2 s _1 (L x = 
6.0 xlO 26 

-1-12 



), the exposure of night 2 corresponds to 3.2x 

,-2„-1/t _ £ a s. i n26 „,.„ „-l\ 



erg s 

10~ 12 erg cirT 2 s -1 (L x = 6.4 x 10 26 erg s -1 ), and the data from 



night 3 give 6.3 x 10~ 12 erg cm" 2 s" 1 (L x = 1.2 x 10 2/ ergs" 1 )- 
The higher value for the third night is caused by the strong flare. 

During all three nights, the X-ray flux in the quiescent time 
intervals is at a level of 2.6 x 10~ 12 erg cm" 2 s (L x = 4.9 x 
10 26 erg s" 1 ). This is comparable to the lower boundary of the 
vari able quiescent X-ra y luminosity 4- 16 x 10 26 erg s _I found 



bv lHaisch etai1dl990h . 



We can calculate the activity level of Proxima Centauri us- 
ing the activity indicator logLx I 'Lboi- With infrared H and K 
band magnitudes of mg = 4.835 and — 4.384 and using 
the bolometric corrections given bv lReid et al. (2001), Proxima 



Centauri's bolometric luminosity is 6x 10 30 erg s ', i.e., slightly 



less th an the Lboi — 6.7 x 10 ergs found by Fro gel et al 



dl972l) . With an activity indicator of logLx /^w = -4 in non- 



flaring time intervals, Proxima Centauri is a moderately active 
star. 

3.2. The large flare 

The X-ray light curve from night 3 (Fig. [3) shows a large 
flare with a peak X-ray flux of 3.8 x 10~ n ergcirT 2 s~' (L x = 
7.2 x 10 27 erg s~'), which is also evident in the OM and optical 
data. The optical UVES exposure-meter and OM light curves 
resemble those of a typical solar impulsive flare. The X-ray 
data and light curves of chromospheric emission lines show 
a more complex behaviour with two broad secondary events 
peaking at about 6:45 and 9:10 UT in the X-ray light curve. In 
the Hor-emission, we find an even more complicated behaviour 
with several sub-peaks in the decay phase, which are not seen in 
any of the other light curves but are roughly associated with the 
two secondary peaks seen in X-rays. The two secondary peaks 
can also be identified in the other Balmer lines and the Can H 
and K lines despite their lower time resolution. The secondary 
peaks are found in neither the optical UVES light curve, nor 
the OM light c urve. A similar e vent was described for Proxima 
Centauri bv lGiidel etail d2002h . 

The flare rise to the maximum flux took place on a 
timescale of ~500 sec. With a 10s binning for X-ray and op- 
tical data, we found that the optical peak as seen in the OM 
light curve precedes the X-ray peak by about two minutes. This 
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Fig. 3. As Fig. Q~| but for night 3. The time intervals for the 
UVES blue arm spectra are also indicated for further reference. 
The dashed line in the OM panel marks the saturation limit that 
was reached in the flare peak. The figure legend is the same as 
in Fig.[T] 
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Fig. 4. PN spectra from night 1 (bottom/red), night 2 (mid- 
dle/blue), and night 3 (top/black). 



can be explained in terms of the Neupert effect (Neupert 1968) 
known to occur in solar and a variety of stellar flares, where 
the time integral of emission due to particle acceleration, such 
as radio emission or blue continuum emission, resembles the 
rise of the flare light curve in soft X-rays. The explanation of 
this behaviour is that accelerated particles hit the dense chro- 




time (h since MJD 54904) 



Fig. 5. Illustration of the Neupert effect during the large flare of 
Proxima Centauri on night 3. Depicted are the combined EPIC 
X-ray light curve in red, its time derivative (smoothed by five 
bins) in black, and the optical UVES light curve in green. 



mosphere and that chromospheric material is heated and evap- 
orates into the corona, where it accumulates and provides the 
emission measure of the subsequent soft X-ray emission. 

This effect is also clearly seen in our data, because the time 
derivative of the combined EPIC (MOS1, MOS2 and PN) X- 
ray light curve during the flare rise (see Fig. [5]l matches the 
shape of the optical light curve; we note that we cannot use 
the OM light curve for this analysis, since the OM data are 
saturated during the flare peak, but we use the light curve that 
is extracted from the UVES's blue-arm photometer. 



3.3. Average magnetic field 

The average magnetic field Bf of Proxima Centauri is mea- 
sured using ou r UVES red arm data by employing a method 
introduced by iReiners & Basril (12006b utilising an absorp- 
tion band of molecular FeH close to 1 fim. We provide a 
brief overview of the m ethod here and refer to the paper by 
Reiners & Basril d2006l) for a more detailed description of the 



procedure. The FeH band contains a large number of isolated 
lines, some of which are highly sensitive to the Zeeman ef- 
fect, while others are not. Since Lande-g factors are unavail- 
able, a direct determination of the Zeeman effect is difficult, 
although there have been so me attempts to mod el FeH molec- 
ular lines theoretically (see lShulvak et al. We measure 
the magnetic field by comparing observed FeH lines to spec- 
tra of stars with known magnetic fields. These template stars 
are GJ 1002 (M5.5, no measured magnetic field) and Gl 873 
(M3.5, Bf ~ 3.9 kG. IJohns-KruU & ValentifeoOOh . Before the 
spectra can be compared, the template spectra have to be scaled 
to match the strengths of the FeH absorption to the strength 
found in our spectra of Proxima Centauri. They also have to 
be adjusted to match in terms of rotational velocity. The ob- 
served spectrum is then modelled as a linear combination of 
the template spectra to determine Bf, thus it is assumed that 
the Zeeman broadening is linear in Bf and that the magnetic 
field strength is distributed similarly in all stars. 
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Fig. 6. Result of the fit for UVES red arm spectrum 5 from 
night 1 with Bf = 400 G ± HOG. Black: slope corrected 
Proxima Centauri data; red: scaled and broadened spectrum of 
magnetic reference Gl 873; blue: scaled and broadened spec- 
trum of non-magnetic reference GJ 1002; green: best fit; Y: 
magnetically sensitive lines; N: magnetically insensitive lines; 
vertical dashed lines border the wavelength regions used for 
the fit. 



We measured the average magnetic field Bf for all three 
nights by using a ^-minimisation to find the best-fit relation 
for the interpolation between the template stars and the UVES 
red arm spectra of Proxima Centauri. The fit was calculated 
in two wavelength regions, 9947.3 A-9961.5 A and 9970.95 A- 
9981.0 A, containing 14 absorption lines of FeH with 7 being 
sensitive and 7 being insensitive to the magnetic field, respec- 
tively. A typical example of our fits is shown in Fig. [6] We com- 
pute a formal lcr uncertainty in Bf by searching the range of 
Bf for which x 1 < X 2 m j„ + 1 while varying the four other param- 
eters. This results in an average uncertainty of ABf = 110 G. 
We note that it is only the statistical error that we use to com- 
pare individual differential measurements because the magneti- 
cally insensitive lines remain constant; thus they provide an ac- 
curate calibration particularly when intercomparing exposures. 
In addition to statistical errors, there are systematic errors and 
we estimate the true uncertainty in the absolute measurement 
of the magnetic field Bf to be of the order of a few hundred 
gauss. Hence, to accurately assess any trend in the data the sta- 
tistical error has to be considered, while the systematic errors 
may offset the whole graph. 

There is a relatively strong correlation between the seeing 
conditions and the derived values of Bf, especially in the sec- 
ond and third night with Pearson's product-moment correlation 
coefficients r of r^i - 0.84 and rm - 0.64. This was elim- 
inated by dividing the Bf measurements by a polynomial de- 
scribing the general trend of the seeing taken from the ESO 
archive. The resulting timing behaviour of Bf is shown in Fig. 
[7] The field h as a mean level o f Bf = 250 G with a statistical 
error of 60 G. Reiners & Basril ( 2008 ) also measured the aver- 
age magnetic field of Proxima Centauri with the same method 
using data taken in April 2004 and found magnetic field values 
in the range 450 G < Bf < 750 G. This change of -350 G over 



Fig. 7. Timing behaviour of the average magnetic field Bf for 
the three observed nights. The onset of the flare in the third 
night is marked by a vertical dashed line. 



five years might well be real, but is close to the limits of our 
uncertainty. Hence, it is unclear whether this difference is real 
or a consequence of the measuring inaccuracy. 

As mentioned above, we can use the statistical error in the 
^-minimisation to compare individual measurements, because 
we are only interested in a differential comparison and not in 
the absolute values of Bf. For the first and second nights, the 
magnetic field is consistent with a constant value within the er- 
ror bars. For the third night, there is a dip in the magnetic field 
with its minimum roughly at the time of the flare onset (Fig.[7j. 
In the context of the reconnection model for flares, a weak mag- 
netic field at flare onset is what one would expect. However, we 
note that the measured fields refer to the photosphere, while the 
field annihilation is assumed to take place in the corona. In ad- 
dition, the slow decay before the flare event is surprising for a 
catastrophic event such as magnetic reconnection, while after 
the flare the magnetic field recovers relatively quickly. In sum- 
mary, it remains unclear whether there is a physical association 
of the magnetic field changes with the flare. 

4. Coronal properties of Proxima Centauri 

4.1. X-ray spectral analysis 

Our X-ray observations of Proxima Centauri can be analysed 
to help us derive several time-variable properties of its corona 
such as temperatures, abundances, and densities, as we investi- 
gate in this section. 

4.1 .1 . Spectral fits and elemental abundances 

To visualise the spectral changes of Proxima Centauri associ- 
ated with the large flare, we show the time evolution of the RGS 
spectra during night 3 in Fig. [8] The mean spectrum is shown 
in the uppermost part of the diagram, while the time evolution 
in the form of the light curve is depicted vertically at the right 
border. The mean spectrum resembles a typical RGS spectrum 
of M-dwarfs with O vm Ly a being the strongest line and other 
prominent lines originating from Ne ix, Ne x, Fe xvn, O vn, and 
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C vi. As expected, the diagram shows a brightening in the spec- 
tral lines with the flare onset. The higher plasma temperatures 
during the flare manifest themselves as a relative brightening 
of lines formed at higher coronal temperatures than the cooler 
spectral lines, for example O vm vs. the O vn triplet. However, 
pronounced changes in the density-sensitive ratio of the O vn 
forbidden to intercombination line cannot be identified in this 
time-resolved plot for the flare and the pre-flare period; they 
are only evident in time-integrated measurements for the flar- 
ing and quiescent states (see section 14.1.31 ). 

There is particular interest in studying elemental abun- 
dances during a flare. Since flares transport chromospheric ma- 
terial into the corona by chromospheric evaporation, coronal 
abundance patterns can temporarily change during flares. We 
determin ed the temperatures and abundances relative to solar 
values ( Grevesse & Sauvalll998l) with an iterative procedure 
of global Xspec fits to EPIC and RGS spectra with VAPEC 
plasma models. In these fits, the abundances and emission mea- 
sure are inherently interdependent, hence we made our infer- 
ences from relative changes in the fit parameters for the three 
observations. The errors we give for the fitted parameters were 
calculated with Xspec's error command. The program Xspec 
calculates the confidence intervals for the desired parameter 
by fixing this parameter at a specific value and varying the 
other free parameters until the best fit is reached. New fixed 
values for the parameter of interest are chosen until the error 
margins that were requested in the error command are cov- 
ered. The errors in these plasma fits are correlated; in particu- 
lar for the emission measure and elemental abundances, there is 
a strong interdependence (see below). Thus, in the calculation 
of Xspec's error margins for an individual parameter, all other 
non-frozen parameters are also allowed to vary, yielding effec- 
tively larger error margins than if one (incorrectly) assumes the 
errors to be uncorrelated. 

When fitting the EPIC spectra, we adopted the following 
method to obtain the abundances: the abundances of oxygen, 
neon, iron, magnesium, and silicon were allowed to vary freely 
and independently, but were constrained to be the same among 
all the VAPEC temperature components. The abundances of all 
other elements were frozen to the solar value because they can- 
not be well-constrained with EPIC spectra since they do not 
have prominent spectral lines at energies where EPIC has a 
large effective area. When fitting EPIC spectra together with 
RGS spectra, carbon and nitrogen abundances also were al- 
lowed to vary. To determine the number of VAPEC components 
with different temperatures, we added new components one by 
one until the addition of the next component did not improve 
the fit significantly. 

In the following sections we used the RGS spectra with 
their much higher wavelength resolution for spectral fitting 
simultaneously to the MOS spectra. We combined the spec- 
tra from the quiescent states of all three nights (for a defini- 
tion of quiescence, see section [3~TT i and compared them to the 
spectra obtained from the large flare during the third night. 
To compare these two states, we defined a fixed temperature 
grid consisting of the values 0.14, 0.4, and 1.0 keV (1.6, 4.6, 
and 11.5 MK). The resulting best-fit^ 2 values of the three- 
temperature VAPEC models along with average individual ele- 



Table 2. Three-temperature fit to the X-ray spectra extracted 
from MOS1, MOS2, RGS1, and RGS2 data allowing individ- 
ual elemental abundance and emission measure to vary on a 
fixed temperature grid. Errors given are lcr errors. 



Parameters 



quiescence large flare 



T\ (keV) 
EMi (10 50 cm- 3 ) 

T 2 (keV) 
EM 2 (10 50 cm- 3 ) 

T 3 (keV) 
EM 3 (10 50 cm- 3 ) 
C 
N 
O 
Nc 
Mg 
Si 
Fe 
red.x 2 
D.O.F. 
logL x (0.2-10 keV) 



0.14 



16 +002 

"• -0.02 



17 +007 

"• 1 ' -0.07 



0.4 



35 +0 - 02 

"■ J -0.02 



1.0 



i 47+0.08 

^ -0.07 



OOi+o.oi 
"• wl -o.oi 

57+ - 04 

"""-0.03 
n 7 q+0.09 

"• ' y -o.m 
45 +0 - 03 

"•^ -0.02 

63 +0 05 

"■"-'-0.03 

30 +0 04 

""'"-0.03 

32 +0 - 07 

"• J -0.07 

32 +0 - 02 

" ,J -0.02 

1.96 
784 
26.70 



1 10 +003 

1 - 1 "-0.03 

66 +0 18 

"•""-0.12 

73 +0 - 21 

—0. 19 

45 +0 03 

" -0.03 

56 +0 04 

""'"-0.04 

54 +0 05 

"• J -0.05 

"•^ 7 -0.05 
n 44+0.02 
"•^ -0.02 

1.877 
561 
27.50 



mental abundances for the quiescent state and the flaring state 
are listed in Table [2] The quiescent state is characterised by 
dominant plasma components at 0.14 and 0.4 keV. During the 
flare, there was a pronounced enhancement of the emission 
measure at 0.4 and 1 .0 ke V. The emission measure at low tem- 
peratures is more weakly constrained than for the quiescent 
state; within the errors, the low-temperature emission measure 
is consistent in both fits. 



4.1.2. The FIP effect 

Similar to the Sun, inactive or low-activity-level stars show a 
normal first ionisation potential (FIP) effect where elements 
with a low FIP are enhanced in the corona compared to ele- 
ments with a high FIP. A reversed pattern - the inverse FIP ef- 
fect - with enhanced high-FIP elements and depleted low-FIP 
elements is frequently found in stars of hig her activity (e.g. 

feooih . 



Brinkman et al l d2001l) , lAudard et al.l (I2003I) ). A FIP effect is 
common among stars with low activity levels (logLx/Lh i < 
-4). Owing to the activity level shown by Proxima Centauri, 
one would expect to observe the inverse FlP-effect. In the mea- 
surements for the quiescent state, the inverse FlP-effect might 
be observed but owing to our measurement errors we cannot 
make a strong claim here, which is illustrated in Fig. [9] The 
actual abundance values can be found in Table [2] During the 
flare, the abundances are slightly higher on average, and the 
iron and silicon abundances rise in relation to the high-FIP ele- 
ments. This is in line with the picture that fresh chromospheric 
material is evaporated into the corona during the flare, which 
causes a change in elemental abundances. In previous X-ray 
observations of Pr oxima Centauri, a similar abundance distri- 
bution was found ( Giide l et al ] l2004h . 

To investigate in greater detail the abundance variations, we 
created spectra by dividing the flare data into several time inter- 
vals with the first four lasting 120s, followed by three intervals 
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10 15 20 25 50 55 

X (A) 

Fig. 8. Time evolution of RGS spectra of Proxima Centauri. Prominent emission lines are labelled in the mean spectrum shown 
in the upper part of the graph. In the right panel, the associated light curve can be found. 



FIP[eV] 

Fig. 9. Normalised elemental abundances to oxy- 
g en abundance relative to solar photospheric values 
((Gr evesse & Sauval 1998)) as a function of the first ioni- 



sation potential (FIP) during quiescence (red) and flare (black). 
Dashed line indicates the solar photospheric abundance. 



of 180s, 240s, and 300s and the remaining data being divided 
into spectral intervals of 600s each. The first spectrum covers 
the flare rise and the following time intervals cover the differ- 
ent phases of the decay. These spectra are fitted with combina- 
tions of APEC models, using two temperature components and 
assuming fixed elemental abundances for the flare plasma but 
leaving the iron abundance as a free parameter. To obtain the 
plasma properties, we use the total emission measure EM, i.e. 



the sum of the emission measures of each temperature compo- 
nent 



EM = YjEMi 



and the flare temperature T is defined as an emission-measure 
weighted sum of the temperatures from each flare component 

_ v TjX EM, 



IS 



EM 



We infer that the iron abundance increases from a level of 
0.30 ± 0.03 directly before the onset of the flare to a maximum 
value of 0.58^ [j, when the flare peak is reached (see Table 
0. As oxygen is another element that produces very strong 
emission lines in our X-ray spectra, we also investigated pos- 
sible oxygen abundance changes during the flare in the same 
fashion, but found no evidence of a similar timing behaviour. 
This indicates that fresh plasma material is evaporated from 
the photosphere and chromosphere that has a different compo- 
sition with a higher iron abundance than the quiescent coro- 
nal plasma. This scenario fits very well with the measured Ha 
asymmetries that indicate movements from the chromosphere 
in the direction of the corona (see section l5l2"V As can be seen 
from line shift measurements in Table|9] material is evaporated 
with a velocity of 44.3 km s _1 in the spectrum lasting from 6:12 
to 6:14 UT. If one makes a conservative estimate that the ma- 
terial is accelerated only at the end of the exposure time, it can 
travel for about 90 seconds or about 4000 km before the mea- 
surement no. 1 of Table [3] starts (6:16:00 UT), which would 
place the material within the corona. For the second measure- 
ment of Ha blue-shifts, the situation is not as clear because 
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Table 4. X-ray counts in the Ovu triplet during quiescence 
(total integration time of 49 ks) and during the flare (total inte- 
gration time of 7 ks). 



line 


counts (quiescence) 


counts (flare) 


resonance (21.6 A) 


225 ± 17 


96 ± 12 


intercombination (21.8 A) 


60 ± 11 


39 ±8 


forbidden (22.1 A) 


130 ± 13 


48 ±8 



the blue-shifts persist for a longer time overlapping with mea- 
surement no. 3 (starting at 6:20:00 UT) of higher coronal iron 
abundance. If most of the material is evaporated at the begin- 
ning of the blue-shift, the material moves for about 210 seconds 
or about 2200 km before the X-ray measurement no. 3 starts, 
again putting the material into the corona. 



4.1.3. Coronal densities 

Electron densities can be inferred from line triplets of He-like 
ions. The RGS energy range contains density-sensitive He-like 
triplets of N vi, O vn, Ne ix, Mg xi, and Si xm. These He-like 
triplets show in increasing order of wavelength a resonance (r), 
intercombination (i), and a forbidden line (f). If the electron 
collision rate is sufficiently high, ions in the upper level of the 
forbidden transition, ls2s 3 S y, do not return to the ground level, 
Is 2 'So, but the ions are instead collisionally excited to the up- 
per level of the intercombination transitions, ls2p 3 P^2, from 
where they decay radiatively to the ground s tate. Thus, the ra- 
tio of the f to i fluxes is sensitive to density ( Gabriel & Jordanl 
(fl%9k 



The He-like triplet of O vn is strong enough in our obser- 
vations to be used to obtain characteristic electron densities in 
the source region. We estimate the coronal plasma densities us- 
ing the density-sensitive ratio of the forbidden to intercombi- 
nation line of the O vn triplet in the quiescent and the flaring 
state. The measured f/i flux ratios are 2.17 + 0.61 in quies- 
cence and 1.23 ± 0.46 during the flare in the third night (see 
Fig- EH Table @}. The errors in these line ratios are rather large 
and are caused by the weak intercombination line. Formally, 
the errors overlap, hence one could argue that the change in 
density is insignificant. However, the ratio (/ + i)/r is prac- 
tically constant during quiescence (0.84 ± 0.17) and the flare 
(0.90 + 0.28). This indicates that physically X-ray photons are 
only shifted from the forbidden to the intercombination line at 
high densities. Hence, the constant (/ + i)/r ratio suggests that 
the change in the f/i ratio is actually density-related and makes 
it less probable that it is a statistical fluctuation. 

To convert the measured f/i ratios into coronal plasma den- 
sities, we approximated the flux ratio by f/i = R a /(1 + n e /N c ), 
where R a is the low density limit and N c is the cr i tical d en- 
sity. We adopted the values from iPradhan & Shulll d 1 98 lh of 
3.95 and 3.1 x 10 10 cm 3 , respectively. We find electron den- 
sities of n e = 2.5 ± 0.7 x 10 10 cirT 3 in quiescence and n e = 
6.9 + 2.6 x 10 10 cm~ 3 during the third-night flaring period. This 
suggests an increase in plasma density during the flare, as one 
expects from the standard flare model. 




A / A 

Fig. 10. Density-sensitive line triplet of O vn (resonance, in- 
tercombination, and forbidden line for increasing wavelength). 
The upper panel shows the best fit to cumulated RGS 1 data dur- 
ing quiescence, while the lower panel shows the RGS 1 triplet 
data during the large flare. 



4.2. Interpretation: loop properties 

Emission measure - temperature (EM-T) diagrams are use- 
ful for estimating physical quantities that are not directly 
observable. Since there is no direct information about the 
morphology of the involved coronal structures, we estimate 
the size of the spatially unresolved stellar coronal flaring 
region from the light curve, time-resolved temperature (T), 
and emission measure (EM) values obtained during the decay 
of the flaring loop. The decay time of an X-ray flare is directly 
related to the length of the flaring loop. However, if there is a 
significant amount of heating present during the decay of the 
flare, it can prolong the decay and lead to an overestimated 
loop length. We assumed that the flare occurs inside a single 
coronal loop with a constant cross-section that is anchored in 
the photosphere. The plasma is confined to the loop, and the 
decay begins after the loop has reached a quasi-steady state, 
where the energy and plasma flow are negligible. According to 
this assumption, the decay time of the X-ray emission roughly 
scales with the plasma cooling time, which in turn scales with 
the length of the loop structure. Therefore, t he longer the decay 
takes, the larger the loop structure becomes. ISerio et alj (1 1 99 If) 
derived a thermodynamic timescale for pure cooli ng of flaring 
plasm a confined to a single flaring coronal loop. Reale et aD 



derived an empirical formula to determine the loop 



length considering the effect of sustained heating during the 
flare decay that uses the trajectory of the flare in the EM-T 
diagram, viz. 
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Table 3. Model parameters for the flare spectra with variable Fe abundances. 



Spect. 


KT, 


EMi 


KT 2 


EM, 


Fe 


red.x 2 d.o.f 


T 


EM 


No. 


[keV] 


[10 5( W 3 ] 


[keV] 


[10 50 'cmr 3 ] 






MK 


[10 50 cm- 3 ] 



1 


75 +005 

"• -0.06 


9 to+0.91 
'" 5S -0 63 


1 90 +0 - 28 

7 -0 26 


4 OH +0-58 
7 " -0.72 


°- 58 -ai2 


1.15 


119 


17.71 + |-« 


7.2s:; : « 


2 


69 +0 07 

"■" -0.09 


3 23 +U2 

J ' -0.85 


1 so+0.19 
1.3»_ 017 


6 39+ - 81 

"• J -0.98 


°- 38 -ao7 


1.10 


148 


14.88+i : 44 


9.62!| 93 


3 


o 6r 003 

"•"-0.03 


9 os+0.53 


1 60 +0 " 

"""-0.12 


4 98 +a4 ° 

t-^o 042 


54+o.n 

-0.09 


1.19 


159 


14 40+ - 82 


7 a +0.94 
'■ 8; '-0.86 


4 


0.63*$ 


3 61 +1 - 47 


1 26 +0 " 27 

"-0 15 


9 pc+1.71 
i,OJ -0.92 


35+ 011 

"• J -'-0 09 


0.95 


93 


10.60 + |- 8 ? 


6 47 +319 

"• -2.80 


5 


0.58!^ 


9 41+0.59 

-0 54 


1 99+O.48 
"' -0.24 


1 Q 1+0.48 
-0.62 


33 +a07 


1.07 


113 


10.05 +2:9 5 


4 09+1.07 
J -1 16 


6 


0.46!™ 


1.805* 


97+ 014 

"•^'-O.B 


1 7 1 +0.38 
" ' -0.75 


24 +0 U 

"■^ -0.07 


0.92 


96 


8.26!j;il 


3 52 +L13 
J-J -1.12 


7 


0.38!»S 


1.16*$ 


1 06 +on 

"""-0.07 


1 67 +0 - 22 


37+ - 20 
"■-"-o.io 


1.20 


106 


9.11!J; 2 ° 


2 84 +0 - 40 

'•° -0.30 


8 


79+° 02 

"• ' -0.03 


9 79+O.34 
'• -0.46 


1 so+0.89 
1.38_ 046 


79+ - 45 

"• ' -0.39 


2i+o.o5 

"•'-0.03 


1.00 


135 


11.28+^ 


Q C9+O.79 
• -0.85 


9 


64 +0 02 

"■" -0.03 


i 7 1 +0.43 
11 ' 1 -0.44 


1 or 012 

""-0.14 


1 09+0.38 
l -° Z -0.35 


26 +0 04 

"•'"-0.03 


1.03 


229 


9 64 +0 82 
y " -0.97 


t co+0.82 
J,JJ -0.79 


10 


62 +0 02 

"•" -0.12 




1 99+O.42 
1,z -0.27 


92 +0 - 28 

"•" -0.32 


30 +0 03 

"• J "-0.03 


1.03 


222 


00+2.27 

8.8»_ 232 


t 70 +0.67 
- 5 -' 8 -0.65 


11 


46 +0 03 


i OO+0.74 
1 - 88 -0.69 


76 +013 

"■ '"-0.06 


2 00 +0 - 72 


27 +0 - 05 

"■"-0.03 


1.00 


217 


7 17 +L22 
'- 1 -0.91 


3 89 +L47 

-1.57 


12 


31+0.03 

" -J -0.04 


i C9+O.29 
"-"-0 32 


72 +0 - 05 

"■ ' -0.04 


1 79+0-34 
■ ' '-0.29 


0.36+™ 7 


0.92 


193 


6 2r - 52 

-0.51 


t 94+O.64 
J " i -0 62 


13 


24 +0 05 

"■' -0.03 


| 9C+O.45 
1-/iJ -0.34 


64 +0 05 

"•" -0.03 


1 70 +0 - 35 

" '"-0.37 


0.38tS 


0.81 


184 


c ca+0.62 
-'• J -'-0.43 


9 qc+o.80 

'■ 7 -'-0.72 


14 


25 +005 

"•'-0.04 


1.04^» 


59+ 06 

-0.03 


1 49+0.38 
1 - H -0.46 


0.38^? 


0.97 


162 


c 9Q+O.73 
J- -0.47 


9 47+0.89 
• -0.82 


15 


22 +0 03 

"'"-0.05 


0-88^ 


57 +0 - 04 

"•-"-0.07 


1 36 +0 - 54 
J -0.29 


34 +0 - 08 

"• J -0.07 


1.14 


142 


5 05 +046 

-'•"-'-0.74 


9 9^+0.83 
• J -0.70 


16 


26 +0 03 

"''"-0.01 


1 67 +0 12 

""'-0.12 


1 7O+0.15 
1 '' 8 -0.09 


40 +0 - 20 

"• H "-0.14 


35 +a21 

"•-' J -0.17 


0.75 


134 


4.28^ 


2.08!JJ; 32 


17 


24 +0 01 
"■'^-o.oi 


1 Q 1+0.10 

i -' 1 -0.11 


0.78!»i 9 


25 +0 - 44 

"•'-0.09 


37+ - 05 

"■-"-0.03 


1.11 


103 


3.82+1 83 


1.57!gS 


18 


24 +0 02 

"■' -0.02 


1 15 +014 

1-1 -0.09 


0.70!™ 


25 +o.ii 


36 +a37 

"■-'"-0.17 


1.02 


91 


t 70+O.83 
-'• '-0.78 


1 40 +0 - 25 

1 - H "-0.24 


19 


31+ 009 

"•■J 1 -0.04 


1 15 +029 


0.68!°| 


21 +0 - 26 

"■'-0.19 


29 +0 07 

"■"-0.07 


0.95 


74 


4 9Q+2.54 
^•"-1.39 


1 ^7+0-56 
"-"-0.50 


20 


22 +0 02 

"•^ -0.02 


79 +0 08 
"•' 7 -o.o9 


69 +014 

"■" -0.09 


15+ 009 

"• 1J -0.06 


20 +0 07 

"■^"-0.11 


1.24 


62 


3 48 +1 01 

-'■^ o -0.61 


95+° i 8 

"• 7 -'-0.16 


21 


13 +003 

"•-0.01 


72 +0 23 

"'"-0.19 


47+ 014 

"■^ ' -0.06 


54 +ai9 

"■ J -0.30 


24 +0 - 05 

"■^*-0.08 


0.98 


60 


t 9C+O.97 
-'•'-'-0.54 


1 26 +0 " 43 

I -'"-0.50 


22 


23 +0 02 

"•' ,: '-0.03 


70 +0 19 

"•'"-0.13 


76 +01 ° 
"•'"-0.10 


45+0.15 

"• HJ -0.13 


32 +ai6 

"•-'■'-0.12 


1.03 


77 


5 14 +a73 

-'• 1 -0.84 


1 15 +a34 

1,1 -'-0.26 


23 


40 +0 08 

"'^"-0.06 


1 05 +021 


091 +0.07 

"■ 71 -o.n 


1 oo +0 - 25 

"""-0.17 


33 +a08 

"■-"-0.05 


1.09 


135 


7 ^7+0.94 
'■-"-1.05 


2 05 +0 - 46 

'■"-'-0.31 


24 


26 +0 07 

"''"-0.04 


Q 76 +0.31 
"•'"-0.23 


n 79+O.O8 
"■"-0.05 


1 99+O.25 
"' -0.28 


n 97+O.O8 
"■ z '-0.06 


0.84 


129 


«r TO+0.90 
6.38_ 059 


1 98 +a57 

" 7o -0.52 


25 


19 +003 


66 +0 13 

"■""-0.14 


0.63!g : « 


0.90«;!? 


n -3 9 +0.09 
"■ J -0.07 


1.14 


105 


5 16 +048 

J - lu -0.44 


1 56 +a32 

i - JO -0.31 


26 


n 9Q+O.O6 
"'' -0.05 


94 +0 22 

-0 22 


0.78!°;1J 


0.38!«j 


32 +a28 

"•-"-0.15 


1.22 


85 


5 09 +1 02 
J -" 7 -i.oo 


1 a a +0.44 
,JJ -0.34 


27 


26 +0 06 

"■'"-0.03 


72 +o:l7 
u./i „ lg 


0.683S 


0.32!°| 8 


32 +018 
"•-"-0.11 


1.02 


69 


4 56 +1 09 

-0.72 


1 04 +036 

1 - KJ -0.29 


28 


23 +0 02 

"'"-0.11 


65 +009 

"■" J -0.16 


o.78:« : ;j 


25 +a21 

"• Z,J -0.08 


25+011 
"•'-'-0.20 


1.14 


66 


4 40+ 1.03 
^■^ 7 -1.53 


91+0-31 
"• 71 -0.25 



J — T ic 



,maxi 



where tic is the decay time derived from the light curve and 
a = 3.7 x 10~ 4 cm^'sK 1 ^ 2 . The observed maximum tempera- 
ture must be corrected to Tq(T - 0, 7) = ^T r> obs (in units of 10 7 
K), where T b s is the maximum best-fit temperature derived 
from spectral fitting to the data. The unit-less correction factor 
is F{() — jf^- + q a , where £ is the slope of EM-T diagram. The 
coefficients t], c a , £ a , and q a depend on the energy response 
of the instrument. According to lReald d2007l) . for XMM/EPIC 
the values are £ = 0.13, 77 = 1.16, c a = 0.51, ( a - 0.35, and 
q a = 1.35. 

We show the EM-T diagram of the flare on Proxima 
Centauri including the two secondary flares in Fig. QT| For 
the part of the flare evolution before the secondary events, the 
slope measured from the logT - 0.51ogEM diagram in Fig.fTTlis 
£ = 1 .07+0.29. Given T peak [MK] = 33.26 +[^3, we determined 
the loop half length to be L = 8.55^ xlO 9 cm. Assuming 
that the volume V oc L 3 , we computed the volume of the flar- 



ing region Vfi are ~ 6.25_jr-4 x 10 cm . We also computed 
the loop-footpoint area V/2L and an area filling-factor of the 
flare of three percent. We also note that in addition to this flare 
decay, there are two other secondary events. From the EM-T 
diagram for these events, the slopes measured from these sec- 
ondary events are similar to the slope obtained for the first flare 




25.2 
EM [(a 



- 3 )] 



Fig. 11. Flare evolution in density- temperature. The green line 
represents the evolution of flare decay (black) and secondary 
events (blue and red) 



event. Thus, the derived loop lengths are of the same order of 
magnitude within the errors. 



4.3. Optical Fe xm emission 

We also investigated the Fexm forbidden coronal line at 
3388. 1 A in our spectra of Proxima Centauri. The stellar spec- 
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trum had previously been searche d for this line before bu t no 
clear conclusion had been drawn (IFuhrmeister et al.ll2004l) . In 
the new UVES data, there is indeed evidence of the forbid- 
den line during and after flare activity for the first and the third 
night. 

To isolate the forbidden coronal line emission at 3388.1 A, 
we averaged at least three normalised spectra for the quies- 
cent state, the flare state, and after-flare state. To emphasise 
the chromospheric and coronal changes between the averaged 
spectra, we subtracted the quiescent state from the flaring state 
and the post-flare state (and for testing purposes also from dif- 
ferent quiescent states of the same night). These difference 
spectra of the different quiescent states usually contain only 
noise as expected. However, the differences between the spec- 
tra of flare states and quiescent states in the first and third nights 
show a weak broad excess at the position of the Fe xm line be- 
ginning at about UT 6:00 for the first night. This line persists 
for the post-flare state. In the third night, the line can also be 
found directly before the large flare starting at UT 4:40, while 
no evidence of Fe xm emission could be found for the second 
night. An example of the line can be found in Fig. [T2J the fit 
parameters obtained with CORA can be found in Table The 
derived fluxes have an error of about a factor of two because 
of the uncertainties in the flux calibration. For the first night, 
we used the spectra no. 3-7 as quiescent state, for the second 
night the spectra no. 4-6, and for the third night the spectra no. 
1-5. The half width cr of the line is much larger than the typical 
measured half width of the chromospheric lines of 0.04 A (see 
Table |6). Interpreting the line width as pure thermal Doppler 
broadening, we expect a half width of 0.04 A, for a typical 
chromospheric temperature of about 10000 K and a half width 
of about 0.17 A for the coronal Fe xm line with a peak forma- 
tion temperature of about 1.6 MK. The measured half widths 
are larger for the forbidden Fe xm line suggesting that there 
are additional line-broadening mechanisms; however, it must 
be kept in mind that the measurement of the line width is sen- 
sitive to the determination of the background, thus the given 
(formal) errors are probably underestimates. 

The order of magnitude of the derived Fe xm flux is con- 
sistent with our estimate derived from Proxima Centauri's 
coronal X-ray spectrum during the flare. The Ovu triplet is 
formed at a similar temperature as Fexm (logT = 6.3 and 
6.2, respectively). We calculate the flux in the O vn triplet dur- 
ing the flare from a CORA fit to the RGS1 spectrum to be 
^Ovu = 4.32 x 10" 13 ergs -1 cm" 2 . The coronal abundance ra- 
tio of iron to oxygen is 0.67 during the flare (see Table and 
from CHIANTI we calculate the ratio of the line emissivities of 
Fe xm to the O vn triplet to be of the order of 0.25 in the tem- 
perature range log T = 6.2 - 6.3. From this, we derive an order- 
of-magnitude estimate of Fp e xm ~ 10~ 15 erg s _1 crrT 2 , which 
is a little bit lower but given the uncertainties compatible with 
the flux computed from the optical spectra in Table [5] 

The observation of the Fe xm line during the larger flare in 
the third night seems to be counter-intuitive at first glance, since 
a flare should heat the corona and Fe xm has its peak ionisation 
equilibrium at T — 1.6 MK. However, during the flare the emis- 
sion measure for the low temperature component at 1.6 MK 
is slightly higher than for quiecent state (or, given the errors, 



4 




3388 
X (A) 



Fig. 12. Example of the fit to the averaged difference spectrum 
of the first night for the after-flare state at the end of the night 
for the Fe xm line. The broad Gaussian is the line best fit ob- 
tained with CORA. 



Table 5. Fit parameters of the Fe xm line. 



averaged 


central wavelength 


half width 


line flux^ 


spectra no. 


[A] 


[A] 


[erg s _1 cm -2 ] 


1st night 


13-17 


3388.02 ± 1.32 


0.25 ± 0.02 


6.610 15 


18-22 


3388.07 ± 1.32 


0.22 ± 0.02 


6.6-10- 15 


23-26 


3388.10 ± 1.32 


0.25 ± 0.02 


7.2-10-' 5 


3rd night 


6-8 


3388.13 ± 1.32 


0.27 ± 0.03 


1.0 -lO" 14 


9-11 


3388.16 ± 1.32 


0.20 ± 0.01 


1.5-10 14 


12-15 


3388.15 ± 1.32 


0.36 ± 0.03 


1.7 -10 14 



" Flux errors are about a factor of two. 



comparable to the emission measure during quiescence) and 
the iron abundance is similary higher, which is consistent with 
the forbidden Fe xm line being seen during the flare. Moreover, 
for another M5.5 s tar, LHS 2076, the line was also observed 
only during a flare (Fuhrmeis ter et alj |2004). 



5. Chromospheric and transition region properties 
of Proxima Centauri 

5.1. Identification of chromospheric emission lines 

For the flare spectrum no. 9 in the blue arm and no. 76 in the 
red arm, respectively, we produced an emission line list con- 
taining 474 chromospheric emission lines, out of which 21 are 
from the red arm. In the blue arm, we analysed the flare spec- 
trum directly, while in the red arm we subtracted the quiescent 
spectrum no. 74. The main part of spectrum no. 9 and spectrum 
no. 7 can be found in Fig. IT3l 

To access the line parameters, we fitted the background, 
central wavelength, half width cr, and line flux as free parame- 
ters using a Gaussian profile for every manually indicated emis- 
sion line in a certain wavelength range. We typically chose 10 
A as the wavelength range size for the fitting process, since the 
background variations for such a short wavelength interval are 
negligible. The line fit parameters including a (possible) iden- 
tification can be found in Table [6] We only show a few rows of 
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Fig. 13. Main part of the UVES blue flare (black) and quiescent (red/grey) spectrum. Several major chromospheric emission 
lines are identified in the plot. 



Table 6. Ten rows of the line catalogue. The whole table is accessible electronically. 



central wavelength 
[A] 


half width 
[A] 


flux^ catalogued wavelength ion 
[erg s~' cirT 2 ] [A] 


multiplet 


id flag 


comment 


3495.35 


0.06 


8.70 


10~ 


15 


3495.285 


Fei 


238 


2 


blend 3495.370 CrII 2 


3495.82 


0.04 


1.25 


10~ 


14 


3495.831 


Mnn 


3 







3496.80 


0.03 


1.01 


10 


14 


3496.814 


Mnn 


3 







3497.09 


0.04 


3.93 


10 


15 


3497.115 


Fei 


78 


2 


cm 


3497.51 


0.03 


7.56 


10 


15 


3497.536 


Mnn 


3 





cm 


3497.83 


0.03 


1.21 


10 


14 


3497.843 


Fei 


6 





cm 


3500.33 


0.04 


8.34 


10~ 


15 


3500.340 


Tin 


6 







3500.87 










3500.852 


Nil 


6 





em, line not fitable 


3503.46 


0.05 


7.53 


io- 


15 


3503.474 


Fen 


4 







3504.89 


0.02 


1.23 


10 


15 


3504.890 


Tin 


88 








" As discussed in section |2~T1 the absolute errors for the given fluxes are about a factor of two. 



the Table [6] in this paper as an example, while the whole table 
is provided only in electronic form. 

The flux measurements given in the line catalogue may be 
affected by rather large errors mainly for the following rea- 
sons. First, the description with a Gaussian may give a poor 
fit quality if the lines have broad wings (and should be fit- 
ted with a Lorentzian/Voigt). Second, the background/quasi- 
continuum may be ill-defined. This is true for emission lines 
in a broader emission line wing, but also for emission cores 
of absorption lines and a variable background. Moreover, nine 
Balmer/Paschen and other lines could not be fitted at all be- 
cause they are too broad or have an ill-defined background level 
(but the lines are clearly there). The absolute flux level is also 
estimated to have rather large errors as discussed in section [2Tl 

For the line i dentification, we generally used the catalogue 
of dMoorelll972l) . A few lines were identified using the NIST 
database^. For the identifications from the Moore catalogue, 
the multiplet is also given in Table [6] The spectra were first 
shifted to the rest wavelength for the identification process 
which should enable us to identify systematic line shifts caused 
by the flare. 



Three lines could not be identified with any known line. 
Another 98 lines have only possible identifications for the fol- 
lowing reasons: (1) the line was not found in the Moore cata- 
logue, (2) the wavelength shift to the possible laboratory wave- 
length is large, (3) the line is blended severely with other lines, 
and (4) (most often) the line is the only one out of the mul- 
tiplet. Since we excluded doubtful features from our line list, 
the list cannot be claimed to be complete, especially for weak 
lines. The identified lines belong in total to 32 different ions. 
The number of lines found for the individual ions is shown in 
Table|7] Statistically, all lines are blue-shifted by 0.01 A. 

We found Balmer lines up to number 19 during as well as 
outside the flare. The lines are only stronger during the flare. 
This indicates that the densities during the flare are not high 
enough to cause any Stark broadening of the lines, whi ch would 
lead to the merging of higher order Balmer lines dSvestka 



19720 . This agrees well with the electron densities measured 



3 Available online under 
http://physics.nist.gov/PhysRefData/ASD/index.html 



using the X-ray data (see Subsection l4.1.3l ), which have at most 
doubled during the flare compared to quiescence. 

We compare the line list t o that of the large flare on 
CN Leo ( Fuhrmeist er et al . 2008), which consists of 1 143 chro- 
mospheric emission lines. The two catalogues overlap in the 
wavelength region from 3280 to 3850 A. In this range, we com- 
pared the proportion of lines found in the CN Leo flare to those 
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Table 7. Chromospheric emission line identifications 



ion 


no. lines 1 


ion 


no. lines 1 


Hi 


19(19) 


Mgi 


10(6) 


He i 


14(6) 


Kn 


1(0) 


Li ii 


1(0) 


Cai 


1(0) 


Sen 


7(4) 


Can 


7(7) 


Til 


3(0) 


Ci 


5(4) 


Tin 


58(53) 


Oi 


3(3) 


Vn 


2(0) 


Nei 


1(0) 


Cri 


9(7) 


Nen 


1(0) 


Cm 


30(26) 


Ah 


2(2) 


Mm 


4(3) 


Si i 


2(0) 


Mnn 


9(8) 


Sin 


4(3) 


Fei 


182(148) 


Cln 


1(0) 


Fen 


45(35) 


Cen 


1(0) 


Co i 


5(2) 


Srn 


2(2) 


Nil 


31(28) 


Zrn 


3(2) 


Nin 


5(5) 


Yn 


2(0) 



in brackets we indicate the number of secure identifications 



found in this event. While for Fe i, Ti i, Mn i, Ni i, and Cr i, the 
proportion of lines found in the CN Leo flare compared to the 
lines found in the Proxima Centauri flare is (partially much) 
larger than 2, for the single ionised species of the same ele- 
ments the proportion is 1—1.5. Moreover, nineteen lines in our 
present catalogue are not found in the CN Leo flare. These are 
4 He i lines, 2 Fe i, 6 Fe n, 2 Zr n lines, 2 Co i, 1 V n, 1 Ti i, and 
1 Hi. This seems to indicate that the flare energy is deposited 
very effectively in the higher chromosphere/lower transition re- 
gion, where helium and Balmer lines are formed, and also the 
lines of singly ionised metals should be formed. 

5.2. Line asymmetries 

During the flare, line asymmetries are found in both the blue 
and the red spectra. The Balmer lines, He r, and Ca n K show 
additional flux in the red wing. In the blue arm, we investi- 
gated the Balmer lines Hy, H6, and the Balmer lines at 3889 A 
and 3835 A. We refrained from fitting He because of the heavy 
blending with the Can H line. We also failed to fit the higher 
Balmer lines because of the extreme width of these lines dur- 
ing the flare, the blending with multiple metal lines, and the low 
line amplitudes of the very highest Balmer lines. Unfortunately, 
all He i lines in the blue arm are weak or blended and decay 
very fast, so that in the spectrum no. 10 after the flare peak 
spectrum no. 9, the helium lines have almost vanished. The 
Ca ii K line being an unblended strong metal line shows a clear 
red wing asymmetry. No other metal lines have clear asymme- 
tries. In the red arm, we find an asymmetry for the He i line at 
6678 A in the flare peak spectrum no. 76, while in spectrum no. 
77, the line has already faded too much to display any asym- 
metry. The other three helium lines in the red arm spectrum 
are either too weak or too complicated in their structure (blend 
with another helium line) to show unambiguous asymmetries. 
Hence, the only line for which we could study the asymme- 
tries with higher time resolution of the red arm spectra is Ha. 



This line shows indeed that the asymmetries have a strong time 
dependence that cannot be studied with the time resolution of 
the blue arm spectra. We fitted the line asymmetries with two 
Gaussians: one for the main (narrow) line component and one 
for the (broad) wing component. For all fits, we subtracted the 
quiescent spectrum directly before the flare. The shifts between 
the broad and the narrow Gaussian are listed in Table |8]for var- 
ious lines and for the Ha line with higher time resolution in 
Table [9] Although the formal errors found in the fitting process 
for the line shifts are small, we assume a ten percent error for 
all the line shifts. 

The high time resolution of He clearly shows that the asym- 
metries are a dynamical phenomenon as can be seen in Fig. [14] 
The asymmetries for the first secondary flare can be identified 
in the time map itself (elapsed time is about 1.3 h). In addition, 
it can be seen in the top panel of Fig. [14] where the yellow 
shaded spectrum shows the most prominent asymmetry in the 
wing. The third peak in this light curve also shows broad wings 
in the Her spectra, which are, however, weaker, thus could not 
be fitted. Starting with spectrum no. 149 (08:52:10 UT) and no. 
150 with a blue asymmetry, the asymmetry turns to the red side 
in spectrum no. 152 and 153 (09:00:48 UT). These asymme- 
tries can also be noted in the time map (elapsed time about 3.2 
h in Fig.lT4T>. While in the blue spectra there are only red asym- 
metries, in Ha there are also blue ones. Indeed for the entire 
integration time of the blue spectrum no. 9, there are mostly 
blue-shifts for Ho-; only for red spectrum no. 76 correspond- 
ing to the first few minutes exposure time of spectrum no. 9 
there is a red shift (the spectrum no. 9 corresponds to elapsed 
time from 0.4 to 0.7 h in Fig. [T4l >. Since the amplitude decays 
exponentially, the conditions at the beginning of the spectrum 
should dominate the spectrum. Therefore, blue spectrum no. 9 
and spectrum no. 76 yield similar red shifts. After red spec- 
trum no. 80 (corresponding to an elapsed time of about 0.6 h), 
the asymmetries and large wings disappear - but re-appear in 
red spectrum no. 97 - 102 (corresponding to an elapsed time 
of about 1.1 to 1.3 h in Fi g . IT4T> . This agrees with the red wing 
asymmetries in blue spectrum no. 10, although the asymme- 
try in the red arm is invisible for the whole integration time of 
this spectrum. That the time scales on which the asymmetries 
evolve can be very short can be seen in red spectrum no. 75 
and blue spectrum no. 8. While the red spectrum shows a clear 
blue-shift, in the blue spectrum the lines either cannot be fit- 
ted or show no real line shifts, but do display broad wings (see 
Tables E IS- 

We interpret all these line shifts as manifestations of mas s 
motions in the chromosphere as e. g. lFuhrmeister et alj (12008I) . 
The high velocity blue-shift at the flare onset corresponds to 
evaporation of chromospheric material, while during the decay 
phase raining down as well as rising material can be found. 



5.3. Theoretical modelling of the flare with PHOENIX 

We modelled the flare spectra with theoretical PHOENIX chro- 
mospheric flare spectra. These spect ra were originally com - 
puted to fit the mega-flare on CN Leo dFuhrmeister et al.l2010h . 
The stars CN Leo and Proxima Centauri have similar photo- 
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Fig. 14. The Ha line variations in time with the quiescence spectrum no. 74 subtracted. The top panel shows three individual 
spectra as marked by the horizontal dashed lines in the time map below. On the right hand panel, we show the light curve of the 
integrated flux of the quiescence subtracted spectrum for the shown wavelength region. The time is measured in elapsed time 
since 5:50 UT. 



Table 8. Line shifts of the broad line component for diverse 
lines from the flare-only spectra. 



line 


velocity 


velocity 


velocity 




[kms- 1 ] 


[km s- 1 ] 


[kms- 1 ] 




spec. no. 8(75) 


spec. no. 9(76) 


spec. no. 10(77) 




05:47:14 


06:14:39 


06:45:12 


Hy ( 4340 A) 




4.8 


12.4 


H<5(4101 A) 


3.9 


4.9 


10.2 


Hi 3889 A 




fit ambiguous 


13.1 


Hi 3835 A 


2.1 


fit ambiguous 


9.4 


CanK 




5.1 


21.4 


He i 6678 A 




6.0 





Table 9. Line shifts of the broad line component for the Ha 
line from the flare-only spectra. 



spectrum 


start time 


velocity 


spectrum 


start time 


velocity 


no. 


(UT) 


[kms" 1 ] 


no. 


(UT) 


[kms" 1 ] 


75 


06:12:16 


-44.3 


97 


06:59:50 


8.7 


76 


06:14:29 


5.3 


98 


07:02:00 


13.7 


77 


06:16:37 


-10.5 


99 


07:04:09 


13.7 


78 


06:18:48 


-9.6 


100 


07:06:19 


33.7 


79 


06:20:57 


-5.9 


101 


07:08:28 


32.9 


80 


06:23:06 


1.4 


102 


07:10:37 


3.7 



spheric properties. While for Proxima Centauri, T e s is about 
3100 K, for CN Leo a r eff of 2900 K was used for the un- 



derlying photosphere in the model calculation. Short & Doyle 
dl998l) studied the influence of the underlying photosphere on 



a quiescent chromospheric model and found that a difference 
of 200 K in T e ff leads to uncertainties in the parameters of the 
column mass at the onset of the transition region as well as the 
column mass of the temperature minimum of about 0.3 dex. 
The finest step in our flaring model grid is 0.2 dex. Since we 
use a photosphere with a lower T e g than is realistic, we should 
infer a too high column mass for the onset of the transition re- 
gion and the temperature mi nimum. Details about th e model 
construction can be found in Fuhrmeist er et al. (2010). In this 
paper, we constructed the flaring model as a linear combination 
of a quiescent spectrum observed directly before the flare and 
a flaring model spectrum. 

For the red and the blue arm respectively, we then used 
a global normalisation of the spectra and a grid of different 
filling factors of the flare model spectra to the quiescent ob- 
served spectra to help us find the best fit to the observed flare 
spectra. The filling factors used for the quiescent spectra are 
0.85, 0.90, 0.93, 0.95, 0.96, 0.97, 0.98, 0.99, 0.995, and 0.999. 
The best fit was found using a x 2 analysis again following 
Fuhr meister et al.l (12010 ). The mean best-fit parameters for the 
five best-fit flare models can be found in TablefTUlwith the stan- 
dard deviation of the mean. Examples of the temperature struc- 
ture of the best-fit models for a selection of different spectra 
can be found in Fig. Q3] For a small wavelength range in the 
blue, a comparison between the observed spectrum and best-fit 
model spectrum can be found in Fig. [TBI 

The model parameters of the fitting process from the blue 
and the red arm agree well with each other with both exhibiting 
a shift towards flare models further out in the atmosphere later 
in the decay phase. For red models later than no. 82, no mean- 
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Fig. 15. Example of the temperature structure of the best-fit 
flare models for different spectra: black line: red spectrum no. 
76 (flare onset); dotted line: red spectrum no. 97 (secondary 
flare); dashed line: blue spectrum no. 9 



flare models. For the higher-pressure models, a higher temper- 
ature for the onset of the transition region is also found. 

The only disagreement between model parameters can be 
found for the filling factor. Comparing red spectrum no. 76 to 
blue spectrum no. 9 shows that the blue spectrum has a lower 
filling factor. This may indicate that there is an error in the nor- 
malisation, which was done for the red arm where the influ- 
ence of the flare on the continuum should be small and then 
applied to the blue arm. Another possibility of the different 
filling factors is a height dependence of the filling factor. The 
blue arm with its wealth of metal lines should trace the mid- 
chromosphere, while the red arm with the Hei, Balmer, and 
Paschen lines should trace the upper chromosphere. The filling 
factor found for the red spectra models of about two percent 
agrees well with that found by the EM-T diagnostics from the 
X-ray data, which is three percent (see section l4~2l ). 



6. Discussion 



Table 10. Mean best-fit flare model parameters in time. 



spec. 


T d „„„ 


log 


log 


giad TR 


filling 


no. 


[K] 


cmassj c i, rom 


cmaSSTtnin 




factor 


76 


8400 ± 300 


-2.5 ± 0.2 


-0.3 ± 0.2 


11.0 ±0.1 


2.0 ± 0.5 


77 


8400 ± 200 


-2.5 ±0.1 


-0.3 ± 0.2 


11.1 ±0.3 


2.0 ±0.1 


78 


8300 ± 200 


-2.7 ± 0.3 


-0.4 ± 0.0 


11.1 ±0.3 


2.0 ±0.1 


80 


8200 ± 300 


-2.7 ± 0.3 


-0.4 ± 0.2 


11.1 ±0.3 


2.0 ±0.1 


82 


8300 ± 300 


-2.9 ± 0.2 


-0.9 ± 0.3 


11.0 ±0.4 


4.0 ±0.1 


95 


8000 ± 100 


-3.0 ±0.1 


-1.1 ±0.3 


11.0 ±0.4 


6.0 ± 1.0 


96 


8400 ± 300 


-2.6 ±0.3 


-0.2 ± 0.2 


11.0 ±0.1 


2.0 ± 0.8 


97 


8400 ± 300 


-2.5 ± 0.2 


-0.3 ± 0.2 


11.0± 0.1 


1.5 ±0.4 


98 


8400 ± 300 


-2.7 ± 0.2 


-0.3 ± 0.2 


11.1 ±0.3 


1.5 ±0.4 


99 


8400 ± 300 


-2.7 ± 0.2 


-0.3 ± 0.2 


11.1 ±0.3 


1.5 ±0.4 


101 


8200 ± 300 


-3.0 ±0.1 


-1.1 ±0.5 


10.9 ±0.3 


4.0 ±0.1 


9 


8500 ± 100 


-2.5 ± 0.2 


-0.3 ± 0.2 


11.0 ±0.4 


0.8 ±0.1 


10 


8400 ± 300 


-2.9 ± 0.2 


-0.8 ± 0.3 


11.0 ±0.1 


0.9 ±0.1 




Fig. 16. Exemplary comparison between observed spectrum 
(black) and our best-fit model spectrum (red) for blue spectrum 
no. 9. 



ingful fit could be found with many models exhibiting nearly 
the same^- 2 . Even later during the secondary flare phase, a fit- 
ting of the flaring models also turned out to be possible again, 
starting at red spectrum no. 95. The secondary flare can be seen 
in the models as a re-shift to higher-pressure chromospheric 



The multi-wavelength character of the observations allows us 
to characterise general behaviour and energetics of the flare for 
various wavelength bands, hence different temperature regimes 
and also cross-check these using different instruments. 



6.1. Flare energetics and excess emission 

The flare exhibited a broadband excess emission in the far 
blue of the UVES spectral arm, which can be easily identified 
in Fig. [13] A similar excess was found for larger flares 
and could in these cases be ascribed to blackb ody emission 
dFuhrmeister et al.ll2008t lHawlev & Fisher! 1 19921) . In the flare 
observed here, it was impossible to fit the excess emission with 
a blackbody model. We nevertheless calculated the excess flux, 
noting that the shape of the spectrum suggests that the bulk of 
the broadband emission is not covered by the UVES spectrum. 
We integrated the flare spectrum up to 3600 A and subtracted 
(1) the quiescent flux, (2) the general flux increase seen in the 
blue arm during the flare, and (3) the chromospheric line flux 
in this wavelength interval for the flux excess calculation. We 
also calculated the total of the chromospheric line fluxes in the 
blue arm and in Ha. We compared the Ho- line flux to the value 
one ob tains using the x method introduced bv lWest & Hawlev 
d2008h . which is based only on the equivalent width of the 
line, the bolometric luminosity, and the distance of the star 
(so one needs no absolute flux calibration). Using the x value 
of 0.176 • 10 -4 for an M6 dwarf, L M of 6.7 • 10 30 erg s -1 , and 
a measured peak EW of 3.4, we evaluated the same flux of 



2.0 TO 12 ergs 1 cm 2 as in our flux calibration (or 1.8 TO 

erg s _1 cm -2 



-12 



using Li,oi of 6.0 • 10 30 erg s~'). The various fluxes 
can be found in Table [TT1 together with the resulting energies 
assuming a filling factor for the flare of two percent for the 
time range of the blue UVES spectrum no. 9 from 6: 14 to 6:45 
UT. 
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Table 11. Various fluxes and flare energetics for different wave- 
length regimes. 



Wavelength band 


Instrument 


Flux^ 

[erg s cirT 2 ] 


Energy 
[erg] 


chrom. em. lines 

Ha emission 
blue excess em. 
(3300-3600 A) 

X-rays 
(0.2-10 keV) 


UVES blue arm 
UVES red arm 

UVES blue arm 
PN+RGS XMM 


1.2- 10" 11 
2.0-10 12 

5.3- 10- 10 

2.2- 10 -11 
±0.2-10-" 


8.5- 10 28 

1.6- 10 28 

3.8- 10 30 

1.9- 10 29 



" The errors in the UVES fluxes are about a factor of two (see 
section [2~TV 



6.2. Ha self-absorption 

A comparison to the flux-calibrated s pectru m of Proxima 
Centauri presented by ICincunegui et aL (2007) shows that the 
Ha line amplitude in our spectra is extremely small during the 
whole three days of our observations. The line amplitude dur- 
ing the fla re is comparable t o the line amplitude of the spec- 
trum from ICincunegui et al. (2007), which was taken during a 
quiescent state. This indicates that Proxima Centauri's chromo- 
sphere was generally in a low activity state during the observa- 
tions. The occurring flare enhanced the Ho- line flux only to a 
level that can also be found during a quiescent state (at other 
times). This fits well with the relatively low X-ray flux found 
in our observations (see section [3TTb . 

The Ha line is heavily self-absorbed during our observa- 
tions. This self-absorption is in most cases asymmetric, but the 
line peak at the red side is higher than at the blue; this peak 
asymmetry should not be confused with the wing asymmetry 
described in section l5T2l of this paper. Two examples can easily 
be noted in the 'flare only' spectra shown in the top panel of 
Fig. [14] while for the flare spectrum the two peaks are about 
of the same height and the line centre therefore symmetric, for 
the two spectra from the decay phase the r ed peak is higher 



than th e blue peak. We note that the spectra of Cincunegui et al 



(2007) have lower spectral resolution hiding a possible self- 
reversal. We therefore compared our data to processed HARPS 
data from the ESO archive facility^, which cover the spectrum 
of Proxima Centauri from 3782 to 6912 A, in the time from 
2004 to 2007. These spectra have not yet been published and 
show a wide variation in the H a line amplitude, with the spectra 
from lCincunegui et al. (2007) at the high end of measured Ha 
line amplitudes and our observations at the low end. Like our 
spectra the HARPS spectra exhibit in most cases a self -reversal 
with more flux in the red side peak. 

We checked the red peak asymmetries using our models, 
which also display Ha self-absorption at the line centre 
that is in all cases symmetric. This would be expected for 
one-dimensional (ID) chromospheric models without mass 
motion, where the self-rever sal is a pure non-lo cal thermal 
equilibrium (NLTE) effect. lAllred et al. d2006l) computed 
ID hydro-dynamical models of M dwarf flares. For various 



evolutionary stages of the flare in their simulations, the Ha 
line profiles show the red peak asymmetries also observed in 
our and the HARPS spectra. In the hydro-dynamical models, 
the peak asymmetry is caused by mass motions during the 
simulated flare. It appears reasonable to assume that similar 
mass motions are also common in the quiescent state to explain 
the observed peak asymmetries. 



6.3. Magnetic field 

The dip in the magnetic field coincides in time with the flare 
in the third night. We estimated the energetics as a consistency 
check, in case a physical connection were possible. The level 
of the magnetic field changes from w 300 + 100 G just before 
the flare to as 100 + 100 G at the time the flare starts in the op- 
tical. So, the nominal change in mean magnetic field is 200 G, 
although the errors are very large. As a consistency check how- 
ever, we can calculate that this corresponds to a released mag- 
netic energy density of 

T] = |j = 1.6 X 10 3 erg cirT 3 . 

In section 14.21 we calculated the flaring volume to be 
y flare ~ 6.25 x 10 29 cm 3 , yielding a total released magnetic 
energy of 1.0 x 10 33 erg. We computed the total X-ray energy 
released during the first 6.5 ks of the flare to be 1 .5 x 10 31 erg by 
integrating the flare and quiescent spectral models from Table|2] 
over the energy range 0.2-10keV and calculating the difference 
between the two. This is two orders of magnitude less than our 
estimate for the total released magnetic energy. However, since 
the detailed loop configuration is not known and we simply es- 
timated the flaring volume as V — L? with L being the loop half 
length, the volume itself can only be an order-of-magnitude es- 
timate. Therefore, we interpret the derived numbers as a gen- 
eral confirmation that the emitted X-ray energy during the flare 
is within the energy budget of the theoretically released mag- 
netic energy. The question, of whether the dip in magnetic field 
is really physically connected to the flare must unfortunately 
remain unsolved. 

Other flare stars of similar spectral type do show higher qui- 
escent activity levels. Although magnetic field measurements 
of mid and late-type M dwarfs are still quite rare, there is a time 



series measurement for the M5.5 star CN Leo dReiners et al 
120071) with mean magnetic field of 2.2 kG covering a large flare. 
A dip in Bf at the time of the flare is also seen in these data, al- 
though there are other field variations of similar amplitude that 
are not associated with flares. Thus, it remains unclear whether 
the change in magnetic field is associated with the flare. 



4 Based on data obtained from the ESO Science Archive Facility 
from program 072.C-0488(E). 



Reiners et al. (2007) found a correlated change in the dis- 
tributions of magnetic field Bf and Ha emission in their data 
of CN Leo. We searched for a similar correlation in our data, 
but found no one, which is unsurprising as the magnetic field 
changes are only small. We note that the Ha emission during 
the course of these observations are quite low (see section 16721 
corresponding to a weak magnetic field. Unfortunately, the 
observation of Rei ners & Basril d2008), which found a stronger 
magnetic field for Proxima Centauri (as described in section 
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13.31 1 does not include the Ha line. 



6.4. Line asymmetries 

Line asymmetries during flares have often been found in high 
resolution spectra mostly for Hi and Hei lines. Examples are 
red line asymmetries for a large flare on LHS 2034 (M6.5) 
( Fuhrmeister et al. 20051). blue line asymmetries for a flare 
on AT Mic (M4) (iGunnet al.lll994h . a nd red line asymme 



tries for flares on AD Leo (M3.5Ve) dCrespo-Chacon et al 



2006) (for further examples see als o references there i n). Fo r 
the CN Leo mega-flare described by Fuhrmeister et al.l (12008b . 
blue-wing as well as red-wing asymmetries have also been 
found. Asymmetries were normally ascribed by these au- 
thors to mass motions, which were partially modelled with 
multiple kernels and compared to the chromospheric down- 
ward condensations (CDCs) observed in the Sun in the case 
of the red asymmetries. During the CN Leo flare, asym- 
metries in the wings of Can lines were also found for the 
first time. Surprisingly, we found asymmetries in the Can H 
and K line for this much smaller flare on Proxima Centauri, 



though ICrespo-Chacon et al.l (120061) found that larger flares 
have stronger asymmetries and found no asymmetry in the Ca n 
H and K line in their flare data (their longest flares lasted about 
30 minutes). 

In contrast to the CN Leo mega-flare we also found blue 
wing asymmetries during flare decay and for the first time a 
re-appearance of the asymmetries during the secondary flare 
phase. The vanishing and re-appearance of the broad wings 
combined with the asymmetries of the Ha line coincides with 
the secondary flare process seen in the X-ray light curve (see 
Fig. [3] and Table |5). When the secondary flare started, the 
broad component vanished, but reappeared some time later. 
The reappearance of the broad wings can also be identified in 
the Ha light curve (see Fig. [3} as the peak after the main flare 
peak. The data again illustrates the dynamic behaviour of the 
asymmetry phenomenon, which provides additional evidence 
of mass motions. 



6.5. Theoretical modelling 

The clear secondary flare found by the flare models for the 
chromosphere strengthens our confidence in this approach 
of mixing quiescent observations with flaring model spectra. 
However, the general c aveats of the method presented by 
Fuhrmeister et al. (2010) still apply: the ID models show (for 
simplicity reasons) a linear temperature rise, which can explain 
part of the difficulties in fitting large wavelength regions with 
the same model spectrum equally well. Moreover, the ID 
models are a very rough description of the horizontally as well 
as vertically highly complex temperature pattern (including 
shocks) of the (solar) chromosphere. Therefore, more physical 
chromospheric models should include hydrodynamic simula- 
tions and dynamic ionisation in either ID or 3D. 




07:00 08:00 

Time (UT) 



Fig. 17. Temporal evolution of flare temperature and emission 
measure. 




Time (100s binned) 



Fig. 18. Light curves with 100s binning and hardness ratio of 
flares observed on Proxima Centauri by XMM-Newton on 12 
August 2001 (red) and 14 March 2009 (black), respectively. 



6.6. Secondary flare events 

The nature of this flare with its two secondary events following 
the first peak in the light curve raises the question of what kind 
of magnetic structure we are dealing with. There are two possi- 
bilities: it could be a single loop that experiences two additional 
heat pulses after the initial flare, or an arcade structure where 
several loops light up in sequence. The heat pulse refers to the 
process where the electron beam is accelerated in the coronal 
part of the magnetic flare loop, then propagates along the mag- 
netic field lines down to the chromosphere, where the electrons 
are fully thermalised. A heat pulse re presents the heatin g of 
the atmosphere by the electron beam. iReale et all (2004) in- 
terpreted an earlier long-duration flare light curve of Proxima 
Centauri as being produced by two loops flaring in sequence, 
since the coronal temperature during the second peak was flat- 
topped and could not be adequately explained by a second heat 
pulse in the same loop. In our data however, the temperature 
and emission measure evolution follows the trend seen in the 
light curves closely as can be seen comparing Fig.fTTland Fig. 
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From our analysis of the EM-T diagram for this flare, we 
derive very similar EM-T slopes for all three flare decays, each 
one indicating that a moderate amount of continuous reheating 
is present. The thus derived loop lengths are of the same or- 
der of magnitude. Even if the errors are rather large, there is 
a decreasing trend discernible within these three lengths; the 
lengths derived from the second and third peaks are compara- 
ble, while the length derived from the first peak seems to be 
larger by a factor of roughly two. This, and that the amount 
of reheating during the decays is very similar, suggests that 
we may see a series of similar loops in an arcade lighting up 
in sequence. However, owing to the large errors in the loop 
lengths, the possibility of a single loop with several heat pulses 
cannot be ruled out completely. A comparison of the X-ray 
and hardne s s ratio light curves for the meg a-flare observed by 
Reale et al.1 d2004 and iGiidel et alj d2004 and the flare pre- 
sented here can be found in Fig. [18] The similarity of the light 
curves of the two events (despite the different count rates) is 
quite remarkable. 

7. Summary and conclusions 

We have presented multi-wavelength observations of the dif- 
ferent activity stages of Proxima Centauri with particular em- 
phasis on a long duration flare. We have tried to construct an 
overall picture of the atmosphere from the photosphere to the 
corona. We have determined a low average magnetic field, cor- 
responding to low average X-ray and Ha emission, with the 
latter being heavily self-absorbed. For the flare, we have com- 
piled a chromospheric emission line list and found asymme- 
tries in the broad wings of Balmer, He i, and Can H & K lines, 
which we ascribe to mass motions, which could explain coro- 
nal abundance changes. From the X-ray data, we have deter- 
mined coronal densities, abundances, emission measures, and 
temperatures. The higher emission measure of the low tem- 
perature component and the higher Fe abundance during the 
flare fit well with the detection of a forbidden optical Fexm 
line during this time. The flare light cur ve is similar to the 
mega -flare for the sam e star described by Real e et al. d2004 
and iGudel et alJ |2004) exhibiting two bumps during the de- 
cay phase. The secondary flares are of a similar loop length 
as the main flare, show line asymmetries, and chromospheric 
flare modelling shows that the first one exhibits similar chro- 
mospheric parameters as the main one. All this indicates that 
the events are not independent of each other, but result from the 
same loop or at least an arcade with several arcade loops ignit- 
ing consecutively. This seems to be a typical flare scheme for 
Proxima Cent auri given, that it h s been observed for the sec- 
ond time, and Reale et al. (2004) noted the similarity to solar 
flares. In contrast to Proxima Centauri, CN Leo, another well- 
observed M5.5 but more active dwarf never showed a compa- 
rable flare cascade during a long-duration flare. On the other 
hand, another highly active dM4.5e star, YZ CMi, also exhib- 
ited a series of secondary flares during a white light mega-flare 
as described by Kowalski et al. (2010bja), who also speculated 
that the flare on YZ CMi originated in a complex arcade with 
a sequence of reconnecting loops. If these interpretations are 
correct, mid-type M dwarfs seems to have a flaring loop geom- 



etry that may be similar to the Sun, which will help to place 
constraints on turbulent dynamo theories. 
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